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1. Introduction 

The top quark is currently the heaviest known 
elementary particle. Thus top physics is a very 
promising field with regard to physics beyond 
the Standard Model (SM). For example, its large 
Higgs coupling makes one hope to learn something 
about the Higgs spectrum in particular, or mass 
generating mechanisms in general. Furthermore, 
its large mass is an important premise for decays 
into non-standard particles. Another exceptional 
property of the top quark is its large decay width 
(rBorn ~ 1-56 GeV) as predicted by the SM. Instead 
of undergoing the process of hadronization, the top 
quark is hardly affected by the non-perturbative 
regime of QCD and decays almost exclusively into a 
bottom quark and a W boson by weak interaction. 

2. Known results 

The 0{as) corrections pi to the decay rate T{t ^ 
bW) amount to —9% of the Born result. This is 
comparable to the expected experimental accuracy 
at an NLC which is around 10%. Thus one should 
make sure that the QCD corrections are reliable, in 
the sense that the series in as converges sufficiently 
fast. This is the main motivation for investigating 
the al corrections to this process. The —9% from 
above may be split into a contribution for Mw — 
(—11%) and the effects induced by a non-vanishing 
W mass (+2%). The electroweak corrections B at 
one-loop level are about 2%. 

For the OUxl) corrections there exists a result 
for Mw = g]. It was obtained by performing an 
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expansion in the limit 

(M2 - M^)/M^ « 1 

and taking into account a sufficiently large number 
of terms in the expansion. The result reads 



rt=rBorn(l- 0.09 -0.02), 



(1) 



where the three numbers correspond to the Born, 
0{as), and 0{a'^) terms, respectively. 

The aim of the calculation of [0 was, on the one 
hand, to perform an independent check of the 0{a^) 
terms at M\y = 0, and, on the other hand, to take 
into account effects induced by a non-vanishing W 
mass at this order. 

3. Method 

While in Q vertex diagrams for t -^ bW were 
computed, in B the top quark decay rate was 
obtained via the optical theorem which relates it 
to the imaginary part of the top quark propagator 



Tt ex Im(i;v -t- Ss)| 



=M? 



(2) 



This means that one should calculate S at the point 
q^ = Mf up to 0{GFa1). An example for a diagram 
that contributes to this order is shown in Fig. nl 
The b quark mass can safely be set to zero, and 
for the moment also Mw = will be assumed. 
The analytic evaluation of diagrams like the one 
shown in Fig. n^ is currently neither available for 
general q^, nor Tor the special case of interest, q^ = 
Mf. The only limiting cases that are accessible 
are g = or Mt = 0. However, asymptotic 
expansions provide an efficient tool to obtain 
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Figure 1. Sample diagram contributing to tlie 
top quark propagator at ©(GpaJ). 



approximate results also away from these extreme 
choices. Their application yields series in q^ jM^ (or 
M"^ /q^), with the individual coefficients containing 
the non-analytic structures in terms of logarithms 
of ji^/q^ and ^ jM"^ (/i is the renormalization 
scale). Employing the analyticity properties of the 
approximated function, one obtains the region of 
convergence for the corresponding series. Within 
this region, the full result can be approximated with 
arbitrary accuracy by including sufhciently many 
terms in the expansion. Examples demonstrating 
the quality of such approximations can be found in 

Applying this strategy to top quark decay, in 
a first step one should compute as many terms as 
possible in the expansion around q^ /M"^ . Note that 
one cannot approach the point q^ = M^ from the 
opposite side (i.e. q^ /M^ > 1), because this implies 
top quarks in the final state which is kinematically 
forbidden. The second step is to extrapolate the 
result from the small-q^ region to the point q^ = 
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Figure 2. q dependence of the [4/4] Pade result 
at 0{GFCts) for different values of 5- 



enormous increase in required computer resources. 
Thus one cannot decide upon the choice of ^ a 
posteriori, but has to set it to some definite value 
from the very beginning. On the other hand, it is 
expected that the behavior of the a^ terms w.r.t. the 
choice of ^ is similar to the 0{as) result. Therefore, 
in |4| the gauge parameter dependence was studied 
in some detail at 0{as). This study was mainly 
based on the stability of the Pade results [tti/ti] 
upon variation of m and n for different values of 
^. Another way to find a "reasonable" choice for S, 
is to study the g^-dependence of the Pade results 
close to the physically relevant point q^ = M^. 



At 0{as) this could be done by explicitly Jhe extrapolation to q = M, is expected to work 



resumming the full series in q'^/M^ [^. At 0(al) 
however, a different strategy was pursued in |j] by 
performing a Pade approximation in the variables 
z = q^/M^ and oj ^ {I - ^/T^l)/{1 + ^/T^^). 

Effects of a non- vanishing W mass can be taken 
into account by applying asymptotic expansions 
w.r.t. the relation M^ ^ q^ ^ M^. fn this way 
one obtains a nested series in M^^/Mf and q^ /M"^ . 
The above procedure can then be applied to each 
coefficient of My^/Mf separately. 

One of the questions one is faced with when 
following this approach is gauge (in)dependence. 
The off-shell fermion propagator is not a gauge 
invariant quantity, and only in the limit q^ = M^ 
the QCD gauge parameter ^ formally drops out. 
Due to the fact that one works with a limited 
number of terms here, gauge parameter dependence 
does not vanish exactly but is only expected 
to decrease gradually as soon as a sufficiently 
large number of expansion terms is included. 
Nevertheless, the claim is that by a reasonable 
choice of the gauge parameter the predictive power 
of the result is preserved. At ©(a^) the calculation 
cannot be performed for general ^ due to the 



best if the variation of the approximating function 
close to this point is smoothest. For some values 
of ^, the q^ dependence of the one- loop result near 
q'^/M^ = 1 is shown in Fig. |[ The results of 
were obtained by setting ^ = 0. The validity of this 
choice is further justified by the perfect agreement 
of the approximation to the exact result at 0{as) 
(see Fig. g and the results in the following section). 
Concerning the technical realization of the 
calculation it heavily relies on automatic Feynman 
diagram evaluation with the help of algebraic 
programs |g[. For details we refer to |Q. 

4. Results 

The result for the top decay rate to 0(al) will 
be written in the following way {y = M^^/M^, 
2/0 = (80.4/f75)2): 



r.=ro(^^(")(y) + ^<5«(y)+(^) <5(^)(,)j, 

(3) 
with S''°\yo) = 0.885.... Following the method 
described in Section S, one obtains S^^^{yo) = 
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Figure 3. T(/) means Taylor expansion of / 
around small y = M^ /M"} up to order n (only 
even orders are shown). The solid line represents 
the exact result for 5^^'. 



—2.20(3) which is to be compared with the exact 
number, reading —2.220 . . .. This demonstrates 
the vahdity and the accuracy of the underlying 
approach. 

At 0{al), the result of ||] reads 

5^'^Hy) = -16.7(8) + 5.4(4) y + y^ (11.4(5.0) 



+ 7.3(1) In y)+0(y') 



3n V^o 



15.6(1.1). (4) 



The first number corresponds to the case of 
vanishing W mass. It agrees perfectly with the 
result of 01 which is —16.7(5). The uncertaint 
quoted in^ is slightly larger than the one of 
but is still negligible as compared to the expected 
experimental accuracy at future colliders. The 
series in M'^/Mf is converging very quickly, and 
the total effect of the Mv^-suppressed terms is 
small. The errors in (jj) are added linearly, yielding 
a bigger uncertainty for the sum than for the Mw = 
result. Adding the errors quadratically instead, 
the uncertainty again would be 0.8. 

5. Estimate for b — > uW 

The integration of Eq. (ph over y from to 1 is 
directly related to the decay rate for b ^ ulD a,t 



0(0^) (and with obvious modifications to the two- 
loop QED corrections to T{fj, -^ evv)). However, 
since the y dependence of 5^^' is only known up 
to y^, it is more promising to pull out 5^ \y) and 
to consider the Taylor expansion of 6'^^'^ /S''^^ which 
is expected to vary only slightly in the relevant y- 
range Isl . On the other hand, if more terms in y are 
included, the Taylor series of (5*^*-'/'^*'°^ becomes ill- 
behaved, and one should directly expand (5^'-' . These 
observations are illustrated for i = 1 in Fig. 0. In 
Q both methods were applied, and the results for 
Tib — > ulv) and r(/x — > evD) agree with the ones of 
Q to about 10%. This constitutes a stringent check 
on both calculations. 

Also a direct application of the method 
described in Section to the (four- loop!) diagrams 
corresponding to 6 ^ ulv and ji —t evv was 
performed in fin], and full agreement with H was 
found. 
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